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ABSTRACT 

A standard operating procedure (SOP) is presented for 
calculating the date of freeze-up on the St. Lawrence River at 
llassena, N.Y. The sOP is based on tvo empirical temperature decline 
equations developed for Kingston, Ontario, and llassena, K.Y. , 
respectively. Input data needed to forecast freeze-up consist of the 
forecast Decenber flov rate and heat flux for the St. Lawrence River 
and the vater teiperature at Kingston, Ontario, on the forecast 
initiation date. The material can provide useful learning activity as 
veil as a research tool. (Author/RH) 
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ST* LAWRENCE RIVER FREEZE-UP 
FORECAST EROCEDURE 



R. A. Adsel 

A standard operating procedure (SOP) is presented for 
calculating the date of free?.e-up on the St+ Lawrence River 
at Massena* N.Y. The SOP is based on two empirical tempera- 
ture decline equations developed for Kingston^ Ont+, and 
Massena^ N.Y,, respectively* Input data needed to forecast 
freeze-up consist o£ the forecast December flow rate and 
heat flux for the St+ Lawrence River and the water tempera-- 
ture at Kingston^ Ont*^ on the forecast initiation date* 
Forecasts are made October 1, October 15, November 1, Novem- 
ber 15j December 1^ and December 15. * 



1 , INTRODUCTION 

A forecast teclinique to predict freeze-up on the St. Lawrence River 
at Hassena^ N.Y,, was developed by the Great Lakes Environmental Research 
Laboratory as part of the Demonstration Program to Extend the Navigation 
Season on the Great Lakes and St+ Lawrence Seaway* Figure 1 portrays the 
area of the St. Lawrence River used in the development of the forecast and 
the specific area for which the forecast is made. The forecast technique 
has been in the development and test piiase for the past 3 years* The pur- 
pose of this memorandum is to present a standard operating procedure (SOP) 
to be used in making operational freeze-up forecasts. The reader is 
referred to an article by Adams (l976) for a complete discussion of fore- 
cast technique development. In instances where Adamses forecast technique 
- modified for use in the present SOP, the modifications are noted here. 



2. FORECAST EQUATIONS 

The frtiezti-up forcciist i.s based upon two tempt?ratiire decline equations* 
Equation (1) is as follows: 

V 

" ^0 '^^ 983' 



where Tj is tiie water temperature; nt Massona^ N.Y.^ T^ is the water tempera- 
ture at Kingston^ Ont. , one travel time earlier^ and (j^ is the given average 
river surface heat flux rate between Kingston and Masseaa. All temperatures 
are in ^C. (The travel time^ i.s the time in days needed for a unit 
volume of wjit<?r to travel from Kingston to Mcisigena+) The river is assumed 
to be well mixed and isothermal during the forciL^ast period (October^ Novem- 
ber^ and December). Freeze-up at Mass'^na Is defined as the time when the 
water temperature at Mas.sena is predicted to be 0,3*C, i.c.^ when T^ = 0+3*C 
At that time* it is probable that an ice cover has formed from shore to 
shore in the forecast areiK It f^hotdd be noted that ^ although the ice cover 
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Figure 2* Area of the St* Lmrence River used in the 
development of the forecast* 



may be lost after the initial freeze-up (due to mild weather)^ it Is the 
initial freeze-up of the winter season that is forecast by the SOP. Slitting 
T, in equation (1) equal to 0*3'*C and solving for T^^ the water temperature 
at Kingston one travel time earlier^ results in the following expression; 

By the use of equation (la)^ T^^ the water temperature at Kingston one 
travel time before freeze-up at Massena^ can be calculated if and F are 
known^ Methods of est?.mating and F are given in Section 3 below. 

Equation (2) is then used to estimate the time needed for the water at 
Kingston to cool from the temperature on the date the forecast is initia- 
ted to T J 



ERIC 



2 



AT = (aT - b)t*^. 



(2) 



where AT is the temperature different^ at Kingston from the date the fore- 
cast is made to one travel time prior to freeze-up at Massena (i^e,, T - Tq) ; 
t is the time in days for the water temperature to drop from T to T^; and 
a, b, and c are constants dependent on the forecast start date. Solving 
equation (2) for t and putting it in logarithmic £orm to facilitate com- 
putation results in equation (2a) below: 

log t = l/c [log (T - Tq) - log (aT - bj]. C2a) 

The number of days, H, from the date of the forecast to the forecasted 
date of freeze-up at Massena is then easily calculated by equation (3), 

H = t + F, (3) 

where t is the time in days for the water at Kingston to cool from T to Tq 
and F is the time needed for the water to travel from Kingston to Massena* 



3* FORECAST INPUT 

Input parameters needed for the solution of equations (la), (2a), and 
(3) are T, F, and Q , The value of T, the daily average water temperature 
at Kingston, is obtained from the Kingston water purification plant* It 
should be noted that the location of the water temperature sensor at 
Kingston was changed from 1000 ft off shore in 50 ft of water to 2800 ft 
off shore in 58 ft of water in the summer of 1970* Freeze-up forecasts 
the past 3 years indicate the change in location does not appear to have 
significantly changed forecast accuracy (Appendix B) . 

The value of F, the travel time, is calculated from the mean December 
forecast flow rate of the river obtained from the U*S. Army Corps of Engi- 
neers, Detroit District. The value of is calculated from regression 
equations of Q and T. The values of used to derive the regression 
equations, unlike those in Adamses report (1975), which are based on heat 
budget analysis of climatological data, are based on solution of the fore- 
cast equations working backward from the actual freeze-up date and F and T 
values* 



4* FORECAST STANDARD OPERATING PROCEDURE (SOP) 

4*1 Initial Forecast and Updates 

The October 1 forecast of freeze-up is usually made October 2 as input 
values must be known for the forecast date*. The forecast will be updated 
on October 15, November 1, November 15, December 1, and December 15^ Jf on 
any of these dates input data are not available, the forecast will be made 
on the first date when input data for the forecast date are available* 
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^+2 Forecast Procedure 

1* Call the water purification plant at Kingston to obtain daily 
average water temperature for the date on which the forecast is 
made CT). 

2* Call the U.S* Army Corps of Engineers* Detroit District* to obtain 
the forecast December flow rate for the St* Lawrence River* 

3+ Calculate travel time CF) from the following formula: 

p ^ 2.1303 X 10^ daTs 
December flow rate* 

where the December flow rate is given in cubic feet per second* 

4* Calculate the value of Q^* the December heat flux rate* from one 
of the following equations: 



A. 


October 1 forecast 


*^t 




113.6 - 28. 7T 


(5) 


b/ 


October 15 forecast 






202.6 - 38. 7T 


(6) 




November 1 forecast 


*^t 




120.1 - 42. 


(7) 




November 15 forecast 


*^t 




-102.8 - 1^ 9T 


(8) 




December 1 forecast 


*^t 




-50.7 - 48. 5T 


(9) 




December 15 forecast 


*^t 




38.1 - 97. 6T. 


(10) 



5* Calculate the value of from equation Cla)* 

X F 

where is the water temperature at Kingston one travel time be- 
fore freeze-up occurs at Massena* 

6* Calculate the value of t from equation (2a)* 

log t = l/c flog (T - Tq - log (aT b)]* 

where T is given from step 1 above an^ Tq is calculated from 
step 5+ The values of a^ b* and c will depend on the date the 
forecast is made as follows: 

7 
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To 


be used with: 




a 


b 




c 


A. 


October 1 forecast 


0. 


00935 


0.04920 


1. 


09886 


B. 


October 15 forecast 


0. 


01055 


0.04061 


1. 


09886 


c» 


November 1 forecast 


0. 


01750 


0.09200 


1. 


09886 


D. 


November 15 forecast 


0. 


02727 


0.12952 


1. 


09886 


E. 


December 1 forecast 


0. 


05180 


0.21750 


1. 


09886 


F. 


December 15 forecast 


0. 


01617 


-0.02665 


1. 


09886. 



7. The predicted freeze-up date is then calculated from equation (-3) , 

N = t + F, 

where t is as calculated from step 6, and F is as calculated from step 3- 
This is a modification from Adams (1976) in which a slightly different 
accounting procedure was used, but the results are virtually the same when 
the forecast freeze-up date is converted to a Julian date. To convert 
these data to the Julian date, use the appropriate conversion factor given 
below: 

A. October 1 forecast 274 + N 

B. October 15 forecast 288 + N 
C- November 1 forecast 305 + N 
D. November 15 forecast 319 + N 

December 1 forecast 335 + N 

December 15 forecast 349 + N, 

On leap years add 1 day to the dates predicted, 

4,3 Sample Calculation 

Suppose the following input data (1-e,, 1 and 2 below) were obtained 
for the October 1, 1974, forecast: 

(1) T is 16a''C 

(2) Forecast December flow is 269,000 ft /s 

(3) Solving for F in equation (4), 
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P^2,1301JLi2!days. 7.9 days 
2.69 X 10 

(4) Solving equation (5) for Q^. for October 1, 

Qj.(Oct 1) = 113.6 - 28.7 T = 113.6 - [(28.7) x (16.1)] = -348 

(5) Solving equation (la) for T^, 

X =03-^=03 -<::341jc_L:91 ^ , 
0 "'^ .983 983 ^ 

(6) Solving equation (2a) for log t, 

log t = 1/c [log (T - Tq) - log (aT •- b)]. 
For October 1 a = 0.00935 

b = 0.04920 

c - 1.09886 

Therefore log t = I [log (16.1 - 3.1) - log (0.00935 x 16.1 - 0.0492)] 
1 • 09886 

log t = Y~^S86 ^^^^ ^^^-^^ " (0.101335)1 
log t - 0.910034 ((1.113943) - (-0.994241)] 
log t = 1.9185. 

(7) Solving equation (3), 

N = t + F = 82.9 + 7.9 = 90.8. 

To convert this number to the Julian date^ choose the appropriate equation 
(in this case the October 1 forecast) » 

274 + N ^= 274 + 90.8 = 364.8. 

Therefore the predicted freeze-up date is December 30.8. 

4.4 Interpretation of Predicted Freeze-Up Date 

Whether this freeze-up date of December 30.8 should be rounded to 
December 31 or truncated at December 30 will be left to the discretion of 
the forecastor as no criteria exist for rounding up or down. 
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The most important thing to keep in mind is that the SOP is an analyti- 
cal tool which must be employed with common sense and imagination and in 
conjunction with as many other forecast aids as are available* For example, 
given an October 1 water temperature of 19*0^^0 (an alltime high) and a pro-- 
jected December travel time of 10*0 days, the method would predict an 
abnormally early freeze-up on December 22.9+ This is so primarily because 
the December heat flux rate is directly proportional to the starting tern" 
perature. If, however, the National Weather Service was calling for above 
normal temperatures in their long^-range forecast for the autumn cooling 
period, it would be foolish to adhere blindly to the results of the forecast 
technique alone* 
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APPENDIX A* FLOW-RATE TRAVEL TIME 



The relationship between forecast December flow rates for the St+ 
Lawrence River and the time needed for a mass of water to travel from 
Kingston to Massena is expressed by the following formula: 

^ 2.1303 X 10^ 
flow rate 

where F = the travel time in days and flow rates are given in thousands 
cubic feet per second. The equation was derived from averages of St. 
Lawrence Seaway cross-sectional areas^ volume flow rate^ and distance 
traveled and is as follows* 

P ^ Area (tth distance (£t) ^ ^^^^^^^ da^ ^ ^^^^ 

flow rate (ft /s) 8.64 x 10 s 
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APPENDIX B* EVALUATION 



Evaluations of freeze-up date forecasts for 9 of the past 10 winters 
by use of the SOP were made (table B.2). Note that the lowest standard 
error of estimate (SE) occurred for the November forecasts • The reason 
for this is not apparent. 



Table S.i* T&nperature$ at Kingston C^C) 



Travel 

— KKTT, NOV, DEC, JAN, Tliac 

In Days 

DATE I IS I 15 I tS I F 



196S-66 


IS, 6 


12,2 


10,0 


7,2 


S,6 




3,3 


9,4 




U,^ 


12, ^ 


11,1 


9,4 


7,8 


6,1 


3,3 


9,8 


! 96 7-68 


1S,0 


U,^ 


11,1 


9,^ 


S,6 


3,9 


2,8 


7,6 


1968-69 


17, & 


|S,0 


12, B 




6,7 


3,9 


1,9 


8, A 


1969-70 


16,7 


15,0 


n,i 


10,0 


6,1 


4, A 


1,1 


8,8 


1970-7 i 


17, a 


17,2 


13, :> 


n,i 


7,2 


4,^ 


2,2 


8,2 


1972-73 


16,7 


13,9 


10,6 


8,3 


6,1 


2,8 


l,t 


7,9 


l973-7ii 


17, S 


16,1 


12,0 


8,3 


6,7 


4,7 


2,2 


8,0 


I974-7S 


16, 1 




11,7 


8,9 


6,4 


4,4 


3,3 


7,6 



Table Forecast of Freeze-^Up Dates a$ Made on Various 

Forecast Initiatzon Dates 



FORECAST INITIATION DATES 

Year Oct, J Oi-t, 15 Nov, 1 Nov, IS Dec, 1 Dec, 15 Observed 

Fr*cic-(iP Oatc 



I96S-66 


D«c, 


30, > 


Jan, 


2,0 




4,S 


Jan, 


5, A 


Jan, 


1,9 


Dec, 


27,7 


Jan* 


9 


1966-67 


J,ln , 


2,8 


O^i-r, 


30,8 


Dec, 


29, S 


Doc, 


2S,9 


Dec, 


24,2 


Dec, 


26,6 


Dec, 


26 


1967-68 


Jari, 


3,2 


Dec, 


31,4 


Jdn^ 


1,0 


Dec, 


27,8 


Jan, 


5,4 


De<;, 


31,2 


Jan, 


2 


1968-69 




26,6 




27,0 


Def , 


2S,S 


DtfC, 


27,1 


Dec, 


27.8 


Due, 


29,7 


Dec, 


2S 


1969-70 


De<, 


Z8,5 


De< , 


26,6 


D^'C, 


30,6 


I>ec* 


2S,4 


iiec. 


30,1 


Dec, 


29,0 


Dec, 


28 


1970-71 


Dec, 


26,8 


0<»^', 


23,fi 




24, S 


DecL 


23,9 


DbC, 


26,5 


Dec. 


30,2 


Dec. 


24 


1972-7 ^ 


Dec, 


:9,4 


Dec* 


29,7 


J<i:i, 


3,0 


D*c. 


31,6 


Dec* 


3K4 


Dec, 


29,9 


D«c, 


30 


1973-74 


Dec, 


27,6 


Dec. 


25, S 


Dec, 


28, S 


Dec, 


31.5 


Dec, 


28,3 


Dec, 


30,6 


Dec, 


30 


1974-75 


Dec, 


n,i 


Dec* 


31,4 


Dec, 


29.7 


Dwc, 


29.4 


Dec. 


30.0 


Dec. 


31,5 


Jan, 


1 



SI 



4,60 



3,SS 



2,86 



2,80 



S,23 



4,92 
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